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Because of recent interest in bats as reservoirs of 


emerging diseases, we investigated the presence of viruses 
in Mystacina tuberculata bats in New Zealand. A novel 
alphacoronavirus sequence was detected in guano from 
roosts of M. tuberculata bats in pristine indigenous forest on 


a remote offshore island (Codfish Island). 


H uman settlement in New Zealand is relatively recent 
compared with that of many countries and extends back 
=800 years to when early Polynesian explorers first arrived 
(7). A wide variety of exotic flora and fauna have since been 
introduced, and major changes to the landscape and ecology 
have occurred, particularly once European settlers arrived 
150 years ago. However, large areas of New Zealand that are 
representative of a prehuman state still remain, in particu¬ 
lar offshore, islands such as Whenua hou (Codfish Island), 
which is situated off the southern coast of New Zealand. Be¬ 
fore human settlement, only 3 species of terrestrial mammals 
were present, all of which were bats. Why nonvolant mam¬ 
mals have been absent is unknown, but this absence was 
probably caused by major extinction events and subsequent 
geographic isolation that prevented recolonization (2,3). 

Only 2 bat species remain in New Zealand and both 
are considered vulnerable (www.iucnredlist.org/). The 
long-tailed bat (Chalinolobus tuberculatus) belongs to the 
family Vespertilionidae and is believed to have arrived 
from Australia =1 million years ago (4). The lesser short¬ 
tailed bat (Mystacina tuberculata) is the sole surviving 
member of the family Mystacinidae in New Zealand (5). M. 
tuberculata bats are believed to have diverged from other 
bat species =20 million years ago (6) and have lived in 
geographic isolation from other bat species until the arrival 
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of C. tuberculatus bats (J). M. tuberculata bats also exhibit 
one of the widest feeding ranges of any bat species and 
occupy a niche similar to rats and mice because these bats 
can walk on the ground by using their wings (7). 

Little is known about the microorganisms present in 
these bats. One study found no evidence of pathogenic 
bacteria or lyssaviruses in them, but reported a Sarcocystis 
sp. (5). Given the recent intense interest in bats as a reservoir 
of emerging diseases and the advent of high-throughput 
sequencing as a virus discovery tool, we investigated the 
presence of viruses in M. tuberculata bats. 

The Study 

Four bat guano samples were collected from known 
roost sites on the remote offshore island of Whenua hou 
(Codfish Island) (46°47'S, 167°38'E), which is situated at 
the southern coast of New Zealand. This small island is 
heavily forested and largely unmodified by humans. It has 
special conservation status for the protection of endangered 
species, and public access is not permitted (9). 

Bat guano was held at 4°C during transport (<48 
h), resuspended in 2 mL phosphate-buffered saline, and 
subjected to centrifugation at 6,000 x g for 5 min. RNA 
was extracted from 400 pL of supernatant by using the 
iPrep PureLink Virus Kit (Life Technologies, Carlsbad, 
CA, USA) and eluted into 50 pL reverse transcription 
PCR molecular-grade water (Ambion, Austin, TX, USA). 
Metagenomic sequencing was then conducted for 1 of the 
samples by using an Ulumina MiSeq Instalment (New 
Zealand Genomics Ltd., Massey Genome Service, Massey 
University, Palmerston North, New Zealand) after a series 
of steps involving DNase I treatment, reverse transcription, 
multiple displacement amplification (QIAGEN, Valencia, 
CA, USA), and Illumina TruSeq library preparation (New 
Zealand Genomics, Ltd.). 

A total of 10,749,878 paired-end sequence reads 
of 250 bp were generated and assembled into contigs by 
using Velvet 1.2.07 (10). Assembled contigs were searched 
for viral sequence by comparison to the nonredundant 
nucleotide database in Genbank (downloaded March 2013) 
by using the nucleotide basic local alignment search tool 
(National Center for Biotechnology Information, Bethesda, 
MD, USA). Forty-six contigs showed similarity to known 
germs Alphacoronavirus sequences (Table). (Raw sequence 
data are available on request from the authors.) 

Virus sequence from the genus Alphacoronavirus was 
confirmed in all 4 original guano samples by using a specific 
reverse transcription PCR based on the metagenomic 
data specific for 582 bp of the RNA-dependent RNA 
polymerase (RdRp) gene (GenBank accession nos. 
KF515987-KF515990). The Rdrp sequence was identical 
in all 4 guano samples, and the closest relative was bat 
coronavirus HKU8 (GenBank accession no. DQ249228; 
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Table. Summary of BLASTn output for contigs from bat guano that show identity to alphacoronaviruses, New Zealand* 


Highest scoring BLASTn hit recorded against 

No. contigs matching Range of contig lengths, _ alphacoronaviruses _ 

Alphacoronavirus gene _ alphacoronaviruses _bp_ e-value _ Nucleotide identity (%)t 


Open reading frame lab 

33 

182-1054 

8 x 10~ 93 

581/828 (77) 

Spike protein 

5 

580-1629 

1 x 1CT 19 

362/551 (66) 

Matrix 

4 

251-840 

2 x 10“" 

532/746 (71) 

Nucleocapsid 

4 

536-890 

6 x 10- 6 

79/109 (72) 

*BLASTn, nucleotide basic local alignment search tool 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_ 
fNo. nucleotides identical between query and database sequence, 

TYPE=BlastSearch&LINK_ 

LOC=blasthome). 



79% nt identity, 426/542). Phylogenetic analysis was 
performed for Rdrp (Figure 1) and for the spike protein, as 
derived from metagenomic data (GenBank accession no. 
KF575176) (Figure 2). We propose that these data support 
identification of a new alphacoronavirus, which has been 
designated as Mystacina bat CoV/New Zealand/2013. 


Conclusions 

The discovery of unknown coronaviruses provides 
information for a model of coronavirus evolution (11) 
and contributes to understanding the process of disease 
emergence, as in detection of Middle East respiratory 
syndrome coronavirus (12). The alphacoronavirus 


AY994055 Feline infectious peritonitis virus 
JN856008 Canine coronavirus strain A76 
JQ408981 Feline infectious peritonitis virus strain DF-2 
NC 002306 Feline infectious peritonitis virus 
JQ408980 Feline coronavirus strain DF-2 R3i 
94 AJ271965 Transmissible gastroenteritis virus 
66 DQ811786 TGEV Miller M60 
97 ** DQ811785 TGEV Miller M6 
** DQ811787 PRCV ISU-1 
gg KC175340 Canine coronavirus strain K378 
FJ938051 Feline coronavirus RM 
HM245925 Mink coronavirus strain WD1127 
99 HM245926 Mink coronavirus strain WD1133 
NC 009988-Rh Bat coronavirus HKU2 
2 AF304460 Human coronavirus 229E 

99 JQ410000 Alpaca respiratory coronavirus isolate CA08-1/2008 
NC 010438 Bat coronavirus HKU8 
96 EU420137 Bat coronavirus IB strain AFCD307 

85 EU420138 Bat coronavirus 1A strain AFCD62 
69 JQ989271 Rousettus bat coronavirus HKU10 isolate 175A 
100 JQ989272 Hipposideros bat coronavirus HKU10 isolate TLC1343A 
JQ989268 Hipposideros bat coronavirus HKU10 isolate TLC1310A 
JQ765566 Human coronavirus NL63 strain NL63/DEN/2008/16 
5 JQ765568 Human coronavirus NL63 strain NL63/DEN/2005/193 

JX524171 Human coronavirus NL63 isolate CBJ123 
°° JQ765564 Human coronavirus NL63 strain NL.63/DEN/2009/14 
JQ765570 Human coronavirus NL63 strain NL63/DEN/2005/235 


NC 009657 Scotophilus bat coronavirus 512 
too i NC 003436 Porcine epidemic diarrhea virus 

i JQ023162 Porcine epidemic diarrhea virus strain attenuated DR13 
^ JX088695 Porcine epidemic diarrhea virus strain GD-B 
JX560761 Porcine epidemic diarrhea virus strain SD-M 
5 KF515987 Mystacina bat CoV/New Zealand/1/2013 

KF515988 Mystacina bat CoV/New Zealand/3/2013 
9/ KF515989 Mystacina bat CoVNew Zealand/2/2013 
KF515990 Mystacina bat CoV/New Zealand/4/2013 
EF065513-Ro-batCov Bat coronavirus HKU9-1 
100 AY772062 SARS coronavirus WH20 
9 DQ22305-Bat SARS coronavirus HKU3-1 

EF065509-PiBatCoV Bat coronavirus HKU5-1 
gg DQ415908-rdrp Human coronavirus HKU1 strain Nil genotype A 

gg AY585228-HCov Human coronavirus OC43 strain ATCC VR-759 
57 AF220295-BCoV Bovine coronavirus strain Quebec 
53 EF424623 Giraffe coronavirus US/OH3/2003 

DQ646405-IVB Infectious bronchitis virus isolate TW2575/98 
95 EU111742-SW1 Coronavirus SW1 

54 FJ376619-BuCoV Bulbul coronavirus HKU11 -934 
FJ376621 Thrush coronavirus HKU12-600 
100 FJ376622-MunCoV Munia coronavirus HKU13-3514 
85 JQ065042-PorCov Porcine coronavirus HKU 15 strain HKU 15-44 


y a-CoV 


y h-cov 


> 


y-CoV 


8-CoV 


Figure 1. Phylogenetic tree showing 
genetic relatedness of RNA-dependent 
RNA polymerase amino acid sequences 
for Mystacina sp. bat coronavirus (CoV)/ 
New Zealand/2013 (shown in boldface) 
with those of known coronaviruses. 
Evolutionary history was inferred for 
183 informative amino acid sites by 
using the maximum-likelihood method 
based on the Whilan and Goldman 
model with gamma distribution in MEGA 
5.05 software (www.megasoftware. 
net). Bootstrap values are calculated 
from 1,000 trees (only bootstrap values 
>50% are shown). Scale bar indicates 
nucleotide substitutions per site. TGEV, 
transmissible gastroenteritis CoV; PRCV, 
porcine respiratory CoV; SARS, severe 
acute respiratory syndrome. 
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Figure 2. Phylogenetic tree 
showing genetic relatedness 
of spike protein amino acid 
sequence for Mystacina sp. 
bat coronavirus (CoV)/New 
Zealand/2013 (shown in 
boldface) with those of known 
coronaviruses. Evolutionary 
history was inferred for 492 
informative amino acid sites by 
using the maximum-likelihood 
method based on the Whilan 
and Goldman + F model 
with gamma distribution and 
invariant sites in MEGA 5.05 
software (www.megasoftware. 
net). Bootstrap values are 
calculated from 1,000 trees 
(only bootstrap values >50% 
are shown). Scale bar indicates 
nucleotide substitutions per 
site. TGEV, transmissible 
gastroenteritis CoV; PRCV, 
porcine respiratory CoV; 
SARS, severe acute respiratory 
syndrome. 
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identified in this study from M. tuberculata bat guano 
is unique in respect to the extreme geographic and 
evolutionary isolation of the host bat species, which along 
with C. tuberculatus bats, has been separated from all other 
mammalian species for =1 million years. 

The current estimate for a common ancestor for all 
coronaviruses is 8,100 bce (13). To be consistent with this 
estimate, Mystacina bat coronavirus would need to have 
been introduced to bats on Whenua hou within the past 800 
years since humans first arrived on this island (given that 
the island had no other mammals before this time) (77) or is 
extant to modern alphacoronavirus phylogenetic radiation 
(genesis and expansion). Apart from humans, only 2 
other terrestrial mammals have ever inhabited Whenua 
hou: the brushtail possum (Trichosurus vulpecula ) and 


the Polynesian rat (Rattns exulans), both of which were 
eliminated from the island in the late 1980s (77); neither 
mammal has been reported as a host of alphacoronaviruses. 
Members of the genus Alphacoronavirus infect only 
mammals. Thus, an avian origin for this vims is unlikely. 

An alternative theory of an ancient origin for all 
coronaviruses has recently been proposed that involves 
an alternative evolutionary molecular clock analysis, 
which places the most recent common ancestor many 
millions of years ago (14). The discovery of Mystacina 
bat CoV/New Zealand/2013 vims could lend support 
to such a theory; despite potentially millions of years of 
isolation, it has diverged relatively little from other extant 
alphacoronavimses, as shown by the close relationship of 
the Rdrp and spike protein genes to those of other extant 
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alphacoronavimses (Figures 1, 2). An expanded survey for 
Mystacina bat coronavims in mammals in New Zealand 
and subsequent characterization of viral genomes would 
provide further insights into the origin of coronaviruses. 

No instances of human zoonotic disease from New 
Zealand bat species have been reported. Only a small number 
of conservation staff handle these bats, these staff use 
standard personal protective equipment and work practice. 
Staff are also offered prophylactic rabies vaccination as a 
precautionary measure, even though New Zealand is free 
from rabies. The genus Alphacoronavirus includes several 
human and animal pathogens, but on the basis of phylogenetic 
data in this study, it is not possible to estimate the risk posed 
by Mystacina bat coronavirus to human or animal health. 

New Zealand is not considered a hot spot for emerging 
infectious diseases. This country is free from many human 
and animal diseases, such as rabies and foot-and-mouth 
disease, and infections with human arboviruses because of 
recent colonization by humans and strict biosecurity border 
controls (75). Thus, indigenous wildlife in New Zealand 
has generally been viewed as an almost sterile and unique 
biosphere. Given detection of this coronavirus, more 
thorough characterization ofthe ecology ofviruses and other 
microorganisms in native wildlife should be considered to 
fulfill conservation needs and further safeguard human and 
domestic animal health against cross-species transmission. 
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